Classical and quantum massive cosmology for the open 

FRW universe 



Babak Vakili 1 * and Nima Khosravi 2 ^ 

1 Department of Physics, Chalous Branch, IAU, P.O. Box 46615-397, Chalous, Iran 
Cosmology Group, African Institute for Mathematical Sciences, Muizenberg 7945, Cape Town, South Africa 

o 

5-1 

Q-r Abstract 



< 
OO 



In an open Friedmann-Robertson- Walker (FRW) space background, we study the classical and 
quantum cosmological models in the framework of the recently proposed nonlinear massive gravity 
theory. Although the constraints which are present in this theory prevent it from admitting the 
flat and closed FRW models as its cosmological solutions, for the open FRW universe, it is not the 
case. We have shown that, either in the absence of matter or in the presence of a perfect fluid, 
the classical held equations of such a theory adopt physical solutions for the open FRW model, in 



a* 



> 
in 



which the mass term shows itself as a cosmological constant. These classical solutions consist of 
two distinguishable branches: One is a contacting universe which tends to a future singularity with 
zero size, while another is an expanding universe having a past singularity from which it begins 
its evolution. A classically forbidden region separates these two branches from each other. We 
then employ the familiar canonical quantization procedure in the given cosmological setting to find 
the cosmological wave functions. We use the resulting wave function to investigate the possibility 
of the avoidance of classical singularities due to quantum effects. It is shown that the quantum 
expectation values of the scale factor, although they have either contracting or expanding phases 
like their classical counterparts, are not disconnected from each other. Indeed, the classically 
forbidden region may be replaced by a bouncing period in which the scale factor bounces from the 
contraction to its expansion eras. Using the Bohmian approach of quantum mechanics, we also 
compute the Bohmian trajectory and the quantum potential related to the system, which their 
analysis shows are the direct effects of the mass term on the dynamics of the universe. 
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1 Introduction 

General Relativity (GR) introduced by Einstein began a renaissance in scientific thought which 
changed our viewpoint on the concept of space-time geometry and gravity. The interpretation of 
gravitational force as a modification of geometrical structure of space-time made and makes this force 
distinguishable from other fundamental interactions, although there are arguments which support the 
idea that the other interactions may also have geometrical origin. Because of the unknown behavior 
of gravitational interaction at short distances, this distinction may have some roots in the heart of 
problems with quantum gravity. Therefore, any hope of dealing with such concepts would be in 
vain unless a reliable quantum theory of gravity can be constructed. In the absence of a full theory 
of quantum gravity, it would be then useful to describe its quantum aspects within the context of 
modified theories of gravity. From a field theory point of view, the gravitational force in GR can 
be represented as a field theory in which the space-time metric plays the role of the fields and the 
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particle that is responsible to propagate gravity is named graviton. Then, naturally in comparison 
to the other field theories, one may ask about the different properties of such a particle. The answer 
to this question is deduced by the linearized form of GR and expansion of the space-time metric g^ u , 
around a fixed background geometry r]^ , as g^ u = r]^ u + h^ v , where is the field representation of 
the graviton. Eventually, it is possible to show that the graviton is a massless spin-2 particle. 

Then, since our knowledge about the behavior of gravity at very long distances is also incomplete, 
a question arises: Is it possible to consider a small nonvanishing mass for the graviton, i.e., a massive 
spin-2 particle? In the first attempts to deal with this question, it seemed that adding a mass term 
to the action may be sufficient. This was done by Fierz and Pauli [1]. However, it was shown that 
by considering the number of degrees of freedom, this model suffers from the existence of a ghost 
field, the so-called Boulware-Deser ghost [2], after studying the non-linear terms. This fact made 
massive gravity an abandoned theory for a while. Recently, de Rahm and Gabadadze proposed a 
new scenario in which they have shown that it is possible to have a ghost-free massive gravity even at 
the non-linear level [3]. That was a positive signal in this area, and the early results in this subject 
have been followed by a number of works that address different aspects of massive gravity [4]. As 
in the case of the other modified theories of gravity, it is important to seek cosmological solutions 
in the newly proposed massive theory of gravity. This is done by the authors of Ref . [5] , who show 
that the existence of some constraints prevent the theory from having the nontrivial homogeneous 
and isotropic cosmological solutions. Indeed, what is shown in Ref. [5] is that, beginning with the 
flat FRW ansatz in the context of massive gravity, the corresponding field equations result in nothing 
but the Minkowski metric. However, by re-examination of the conditions, the authors of Refs. [6, 7] 
have shown that, for the open FRW model, this is not the case, and the nonlinear massive gravity 
admits the open FRW as a compatible solution for its field equations. Another progresses to find 
the massive cosmologies lie in the field of the bi-metric theories of gravity; see, for instance, Ref. [8] 
based on the works of Hassan and Rosen [9], in which they show that a bi-metric representation for 
massive gravity exists. 

Our purpose in the present paper is to continue the works of the authors of Ref. [6, 7] in 
greater detail, based on the Hamiltonian formalism of the open FRW cosmology in the framework of 
massive gravity. We obtain the solutions to the vacuum and perfect fluid classical field equations and 
investigate their different aspects, such as the roll of the graviton's mass as a cosmological constant, 
the appearance of singularities, and the late time expansion. We then consider the problem at hand 
in the context of canonical quantum cosmology to see how the classical picture will be modified. 
Our final results show that the singular behavior of the classical cosmology will be replaced by a 
bouncing one when quantum mechanical considerations are taken into account. This means that the 
quantization of the model suggests the existence of a minimal size for the corresponding universe. We 
shall also study the quantum model by the Bohmian approach of quantum mechanics to show how 
the mass term exhibits its direct effects on the evolution of the system. 

The structure of the paper is as follows. In section 2, we briefly present the basic elements of the 
issue of massive gravity and its canonical Hamiltonian for a given open FRW universe. In section 3, 
classical cosmological dynamics is introduced for the vacuum and perfect fluid. Quantization of the 
model is the subject of section 4, and in section 5, the Bohmian approach of quantum mechanics is 
applied to the model. Finally, the conclusions are summarized in section 6. 



2 Preliminary set-up 

In this section, we start by briefly studying the nonlinear massive gravity action presented in Refs. 
[6, 7] for the open FRW model, where the metric is given by 



ds 2 = g^dx^dx" = -N 2 {t)dt 2 + a 2 (t) 



dx 2 + dy 2 + dz 1 



\K\ (xdx + ydy + zdzY 
I + \K\(x 2 + y 2 + z 2 ) 



(1) 
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with N(t) and a{t) being the lapse function and the scale factor, respectively, and K = — 1 denoting 
the curvature index. Here we work in units where c = fi = 16irG = 1. In the massive gravity scenario 
one considers a metric perturbation as [5, 10] 

9tiv = rjiiv + h^u = 7] ab d ll (t) a {x)d 1/ (p b (x) + H^, (2) 

where r] a b = diag(— 1, 1, 1, 1) and (j) a {x) are four scalar fields known as Stiickelberg scalars and are 
introduced to keep the principle of general covariance also in massive general relativity [11]. It is clear 
that the first term in (2) is a representation of the Minkowski space-time in terms of the coordinate 
system {(fp,4> 1 ) and thus the tensor H^ u is responsible for describing the propagation of gravity in 
this space. The action of the model consists of the gravitational part S g and the matter action S m as 

S = S g +S m . (3) 

The matter part of the action is independent of the massive corrections to the gravity part. Also, the 
gravity part can be expressed in terms of the usual Einstein-Hilbert, with an additional correction 
term coming from the massive graviton; that is [5] 

d 4 x, (4) 

in which all of the modifications due to the mass and also the interactions between the tensor fields 
H^ v and g^ u are summarized in the potential U(g,H). By using ghost-free conditions for the theory 
in Ref. [11], we propose the following form for the potential term: [10] 

U{g,H) = -4(£ 2 + a 3 £ 3 + a 4 £ A ), (5) 

where 



R 



^U(g,H) 



■ £ 2 = \ (< /C> 2 - < K? >) , 

< £ 3 = i (< /C> 3 -3 < K >< K 2 > +2 < /C 3 >) , (6) 

k U = ^ (< K. > 4 -6 < K > 2 < K? > +3 < K? > 2 +8 < K >< /C 3 > -6 < /C 4 >) , 
in which the tensor K,^ is defined as 

Kltig, H) = 8$- y/vabd^dvtf, (7) 

and the notations < JC >= g^K,^, < K? >= g 010 g^lC^lCp u ,... are used for the corresponding traces. 
Now, equations (4)-(7) describe the gravitational part of the action for a massive gravity theory. Since 
its explicit form directly depends on the choice of scalar fields (f) a (x), it is appropriate to concentrate 
on this point first. Interesting forms for such fields should involve terms which would describe a 
suitable coordinate transformation on the Minkowski space-time. In a flat FRW background, for 
instance, one may select <fP = f(t) and <p l = x l , as is used in [5]. Here, for the open FRW metric (1), 
we use the following ansatz proposed in [6] 

4> Q = f{t)^l + \K\ Xl x\ ^ = J\K\f(t)x\ (8) 

Upon substitution of these scalar fields and also the definition of the Ricci scalar into the relations (4)- 
(7), we are led to a point-like form for the gravitational Lagrangian in the minisuperspace {N, a, /} 

as 

C g = — ^- - 3\K\Na + m 2 (L 2 + a 3 L 3 + a 4 L 4 ) , (9) 
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where 



L 2 = 3a (a - vTO) (2iVa - a/ - AT^/) , 
L 3 = (a - V^/) 2 (4iVa - 3a/ - iV/TO) , 
L 4 = (a-/^/) 3 (iV-/), 



(10) 



in which an overdot represents differentiation with respect to the time parameter t. It is seen that this 
Lagrangian does not involve N, which means that the momentum conjugate to this variable vanishes. 
In the usual canonical formalism of general relativity, we know this issue as the primary constraint in 
the sense that the variable N is not a dynamical variable but a Lagrange multiplier in the Hamiltonian 
formalism. On the other hand, Lagrangian (9) seems to show an additional constraint related to the 
Stiickelberg scalars whose dynamics are encoded in the function f(t). We see that in spite of the 
common Lagrangians in which the first derivative of the configuration variables are of second order, 
/ appears linearly in the Lagrangian (9). Therefore, by computing the momentum conjugate to /; 



that is, Pf = ^j-, we obtain 



df 



P f = -m 2 (a - y/\K\f) \sa 2 + 3a 3 a(a - y/\K\f) + a 4 (a - J\K\f? 



(11) 



Now, it is clear that this relation is not invertible to obtain f(f,Pf). In such a case, the Lagrangian 
is said to be singular and the relations like (11), which hinder the inversion, are known as primary 
constraints. One may use the method of Lagrange multipliers to analyze the dynamics of the system 
by adding to the Lagrangian all of the primary constraints multiplied by arbitrary functions of time. 
However, to deal with our constrained system, we act differently and proceed as follows. We vary the 
Lagrangian (9) with respect to / to obtain 



KIN 



K\(a 3 + a 4 )r(i) - 2J\K\(1 + 2a 3 + a A )a{t)f{t) + (3 + 3q 3 + a 4 )a 2 (t) 



0. (12) 



The solution a = y/\K\N of this equation is nothing but what we obtain from the variation of the 
usual Einstein-Hilbert Lagrangian with respect to N. Since its counterpart in massive gravity is 



3\K\(a 3 + a 4 ) 2 + m 2 a 2 (t) [2(1 + a 3 + a§ - a 4 ) 3 / 2 - (1 + a 3 )(2 + a 3 + 2a\ - 3a 4 )] 

V3(a 3 + a 4 ) 



N, (13) 



we cannot accept the relation a = \J\K\N as a physical solution. Therefore, the constraint corre- 
sponding to the dynamic of f(t) shows itself in the equation 



\K\(a 3 + a A )f\t) - 2yJ\K\{\ + 2a 3 + a 4 )a(t)f(t) + (3 + 3a 3 + a A )a 2 {t) 
where using the same notation as in [6] , its solutions can be written as 

1 1— 2— -l - . 4- . r 



0, 



(14) 



/(*) 




] + a 4 ± y 1 + a 3 + CK3 — a 4 
a 3 + a 4 



(15) 



As is argued in [6], in the limit where a 3 and q 4 are of the order of a small quantity e, the expression 
of X + goes to infinity while X_ — > 3/2. Because of this limiting behavior, we use the subscript — in 
the following for numerical values of constants with subscript ±. Now we may insert the constraints 
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(15) into the relations (10) to reduce the degrees of freedom of the system and obtain a minimal 
number of dynamical variables. If we do so, we obtain 



L 2 = 3 (1 - X±) 



l 3 = (i - x ± y 



l 4 = (i - x±y 



(2-X±)N 



a 3 , 



(4-X±)N-3 



(16) 



iv-4±=d 



« 3 , 



in terms of which the Lagrangian (9) takes its reduced form with only one physical degree of freedom 
a. The momentum conjugate to a is 



P -2£s. 

oa 



6aa 



Noting that 



dL 2 
da 

one gets 



= 3 



(Xt - l)a 3 



<9L 3 
da 



dL 2 



dL, 



dL 4 



oa oa oa 



(X ± - 1)V, 



\K 



da 



\f\K 



3^3 



(X ± - iya 



Pa = -6— " 



ad C±m 2 



C± = X±{1 -X±) 3 + 3a 3 (l - X±) + a 4 (l - X±f 



(17) 



(18) 



(19) 



Now, the Hamiltonian of the model can be obtained from its standard definition H = aP a — £, with 
result 



H g = NU 9 = N 



1 ( Pa + 9^ a A +3\K\a + c±m 2 a 3 



12a \ " y/\K 
in which we have defined 

c± = {X± - 1) [3(2 - X±) + a 3 (l - X±)(A - X±) + a 4 (l - X ± ) 2 



(20) 



(21) 



We see that the lapse function enters in the Hamiltonian as a Lagrange multiplier as expected. Thus, 
when we vary the Hamiltonian with respect to N, we get H g = 0, which is called the Hamiltonian 
constraint. On a classical level this constraint is equivalent to the Friedmann equation, wherein our 
problem at hand can be easily checked by comparing it with the equation of motion (4.5) in [6]. 
On a quantum level, on the other hand, the operator version of this constraint annihilates the wave 
function of the corresponding universe, leading to the so-called Wheeler-DeWitt equation. 

Now, let us deal with the matter field with which the action of the model is augmented. As we 
have mentioned, the matter part of the action is independent of modifications due to the mass terms. 
Therefore, the matter may come into play in a common way and the total Hamiltonian can be made 
by adding the matter Hamiltonian to the gravitational part of (20). To do this, we consider a perfect 
fluid whose pressure p is linked to its energy density p by the equation of state 



p = up, 



(22) 



where — 1 < oj < 1 is the equation of the stated parameter. According to Schutz's representation for 
the perfect fluid [12], its Hamiltonian can be viewed as (see [13] for details) 



N- 



P 



,3u) ! 



(23) 
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where T is a dynamical variable related to the thermodynamical parameters of the perfect fluid and 
Pt is its conjugate momentum. Finally, we are in a position in which can write the total Hamiltonian 
H = H„ + H m as 



H = NU = N 




+ 3\K\a + c±m l a A + 



Pr_ 

„3u> 



(24) 



The setup for constructing the phase space and writing the Lagrangian and Hamiltonian of the model 
is now complete. In the following section, we shall deal with classical and quantum cosmologies which 
can be extracted from a theory with the previously mentioned Hamiltonian. 



3 Cosmological dynamics: classical point of view 

The classical dynamics are governed by the Hamiltonian equations. To achieve this purpose, we 
divide this section into two parts. We first consider the case in which the matter is absent, i.e., the 
vacuum, and then include the matter. 

3.1 The vacuum classical cosmology 

In this case, we can construct the equations of motion by the Hamiltonian equations with use of the 
Hamiltonian (20). Equivalently, one may directly write the Friedmann equation from the Hamilto- 
nian constraint H = which, as we mentioned previously, reflects the fact that the corresponding 
gravitational theory is a parameterized theory in the sense that its action is invariant under time 
reparameterization. Noting from (19) that 

N f n C±m 2 o\ 



equation (20) gives 

3aa 2 - 3\K\a = c±m 2 a 3 , (26) 

in which we have chosen the gauge N = 1, so that the time parameter t becomes the cosmic time 
r. As is indicated in [6], this equation looks like the Friedmann equation for the open FRW universe 
with an effective cosmological constant A± = c±m? and admits the following solutions 



a ± (r) = A / T sinh ±W-(r-r,) , (27) 




where r* is an integration constant and we have taken A = c_m 2 . For a positive r*, the condition 
a(r) > implies that the expressions of a+(r) and a_ (r) are valid for r > r* and r < — r* respectively, 
such that a±(r*) = 0. It is seen that the evolution of the corresponding universe with the scale factor 
a + (r) begins with a big-bang-like singularity at r = r* and then follows an exponential law expansion 
at late time of cosmic evolution in which the mass term shows itself as a cosmological constant. For 
a universe with the scale factor a_(r), on the other hand, the behavior is opposite. The universe 
decreases its size from large values of scale factor at r = — oo and ends its evolution at r = — r* with 
a zero size. In figure 1 we have plotted theses scale factors for typical values of the parameters. As 
this figure shows, although the behavior of a + (r) (a_(r)) is like a de Sitter (a(r) ~ eV / A^ r ) universe 
at r — > oo (r — > — oo), in spite of the de Sitter, it begins (ends) its evolution with a singularity. 
In summary, what we have shown previously is that in the framework of an open FRW background 
geometry, the vacuum solutions of the massive theory are equivalent to the solutions of the usual GR 
with a cosmological constant. Accordingly, the zero-size singularity of both theories has the same 
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Figure 1: The figures show the evolutionary behavior of the universes based on (27). We have used the numerical 
values A = 1 and r* = 0. 



nature. In this sense we would like to emphasize that the metric (1) with the scale factor (27) is 
indeed a section of the de Sitter hyperboloid 

-T 2 + X 2 + Y 2 + Z 2 + W 2 = 1, (28) 

embedded in a 5-dimensional Minkowski space 

ds 2 = -dT 2 + dX 2 + dY 2 + dZ 2 + dW 2 . (29) 

To see this, one may parameterize the hyperboloid in terms of the spherical coordinates (r, 9, 4>) as 
[14] 



' T = 


V 1 + r 2 sinh 


T, 


X = 


cosh r, 




< Y = 


r sinh r cos <p 


cos 


Z = 


r sinh r cos 4> 


sin 


. W = 


- r sinh r sin (j. 


>, 



(30) 



which, upon substitution into the metric (29), yields the open FRW metric with the scale factor 
a(r) = sinhr. This means that the point a = can be viewed as a coordinate singularity. However, 
we have to note that in the presence of any kind of matter field the point cl{t*) = represents a 
true singularity. Thus, our following analysis to quantize the model is based on the minisuperspace 
coordinate system in terms of which the dynamical representation of the metric, i.e. (1), is written. 

In the next section, we shall see how the previous picture may be modified when one takes into 
account quantum mechanical considerations. 

3.2 Perfect fluid classical cosmology 

Now, we assume that a perfect fluid in its Schutz's representation is coupled with gravity. In this 
case the Hamiltonian (24) describes the dynamics of the system. The equations of motion for T and 
Pt read as 

f = {T,H} = ^j, P T = {P T ,H} = Q. (31) 
A glance at the above equations shows that with choosing the gauge iV = a 3w , we shall have 



N 



,3a> 



(32) 
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which means that variable T may play the role of time in the model. Therefore, the Friedmann 
equation H = can be written in the gauge N = a 3uJ as follows 



3h 2 



3\K\a 6ul + Aa eu}+2 + R 



3w-l 



(33) 



where we take -P7 1 = Po = const, from the second equation of (31). Since it is not possible to find 
the analytical solutions of the above differential equation for any arbitrary oj, we present its solutions 
only in some special cases. 



oj = — i: cosmic string. In this case we obtain 



a(t) 



A 



(t - t f 



P + 3\K\] 1/2 



A 



(34) 



where to is an integration constant. We see that the evolution of the universe based on (34) has 
big-bang- like singularities at t = to ± t* where t* = V^ p "^\ ^jl^ Indeed, the condition a 2 (t) > 
separates two sets of solutions a^(t) and a^ IT \t), each of which is valid for t < to — t* and t > to + i*, 
respectively. For the former, we have a contracting universe which decreases its size according to 
a power law relation and ends its evolution in a singularity at t = to — t*, while for the latter, the 
evolution of the universe begins with a big-bang singularity at t = to + 1* and then follows the power 
law expansion at late time of cosmic evolution. 

One may translate these results in terms of the cosmic time r. Using its relationship with the 
time parameter t in this case, that is, dr = a~ 1 (t)dt, we are led to 



a(r) = I 



a«(r) 



t (")(r) 



_L 

/12A 



_L 

/ I2A 



-(r-ro) _ 3(p Q + 3|^|) e V 3(^0) 



,71 (t-to) _ 3 (p o + 3 \ K \} e -^ (r-ro) 



T < TQ - r*, 



r > T + T*, 



(35) 



where r* 



^1 ln3(Po + 3|-fiT|). Again, it is seen that there is a classically forbidden region To — r* < 
r < to + t*, for which we have no valid classical solutions. For r < To — r*, the universe has 
a exponential decreasing behavior which ends its evolution in a singular point with zero size at 
r = to — t* , while in the region r > To + r* it begins with the big-bang singularity at r = to + t* and 
then grows exponentially forever. 

• oj = —1: cosmological constant. Performing the integration, we get the following implicit 
relation between t and a(t): 



VS(P + A) 2 



6\K\ + (P + AK V 3 I K I + ( p o + A)a 2 = t-t . 



(36) 



In terms of the cosmic time r, it is easy to see that this solution returns to (27), in which the cosmo- 
logical term is replaced by A — > A + cons. This is expected because the solutions (27) were equivalent 
to an open FRW universe with a cosmological constant. Therefore, adding a new cosmological term 
(a perfect fluid with oj = — 1) only makes a shift in the corresponding cosmological constant. 



4 Cosmological dynamics: quantum point of view 

In this section we look for the quantization of the model presented above via the method of canonical 
quantization. As is well known, this procedure is based on the Wheeler-DeWitt equation = 0, 
where % is the operator version of the Hamiltonian constraint and * is the wave function of the 
universe, a function of the 3-geometries and the matter fields. As in the case of the classical cosmology, 
we consider the matter of free and perfect fluid quantum cosmology separately. 
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Before going to the subject, a remark is in order related to the Hamiltonians (20) and (24). The 
term in the round bracket in these Hamiltonians is like the Hamiltonian of a charged particle moving 
in an electromagnetic field. Prom this analogy, one may define the transformation 



n„ 



„ C±m 2 q 
P a + -^—a 3 , 



a, 



(37) 



to simplify the form of the classical Hamiltonian. It is clear that this is a canonical transformation 
both classically and quantum mechanically [15]. Since going back from a new set of variables to the 
old ones in a classical canonical transformation can be made without any ambiguity, applying this 
transformation may not be important for the classical dynamics presented in the previous section. 
In the context of quantum mechanics, on the other hand, the subject is of little difference. The 
transition to the quantum version of the theory is achieved by promoting observables to operators 
which are not necessarily commuting. Thus, by replacing the canonical variables (a, P a ) by their 
operator counterparts (d,P a = —id/da), we obtain the quantum Hamiltonian 



*=—»- 



1x4 2 1 ~ 

+ - = -12 a 



Pa + 



C±m 2 



+ ..., 



(38) 



where ... denotes the terms out of the round bracket in expressions (20) or (24). When calculating the 
square, it should be noted that the operators a and P a do not commute. Although the order of these 
operators does not matter in the classical analysis, quantum mechanically this issue is quite crucial. 
Indeed, this is the operator ordering problem and, unfortunately, there is no well defined principle 
which specifies the order of operators in the passage from classical to quantum theory. There are, 
however, some simple rules which one uses conventionally. If, for instance, we order the products of 
a and P a in n a such that the momentum stands to the right of the scale factor, we obtain 



n„ 2 



Clm 4 
\K\ 



a° + 2 



C±m? 



d 6 P a 



C±m 2 n 

6i — - — -a , 



(39) 



in which we have used the commutation relation [a, P a ] = i. With this expression at hand, there 
is still another factor ordering ambiguity in the terms a~ 1 P a and d 2 P a to construct the quantum 
Hamiltonian (38). As Hawking and Page have shown [16], the choice of different factor ordering will 
not affect semiclassical calculations in quantum cosmology, so for convenience one usually chooses 
a special place for it in the special models. However, in general, the behavior of the wave function 
depends on the chosen factor ordering [17]. In what follows, as one usually does in the minisuperspace 
approximation to the cosmological models, we work in the framework of a special factor ordering in 

a 2 ~ 2 

which, in addition to the expression (39) for YL a , we also use the orderings a~ l P a = P a d~ 1 P a and 
d 2 P a = dP a d to make the Hamiltonian hermitian 1 . 

4.1 The vacuum quantum cosmology 

In this case, with the help of the Hamiltonian (20) and use of the abovementioned choice of ordering, 
the Wheeler-DeWitt equation reads 



da 2 



+ 



da 



36 + 2i— A 



a 2 + 12Aa 4 



a 2 

4 



*(a) = 0. (40) 



1 With the canonical transformation (37) at hand, one may uses the transformed Hamiltonian H — — j^^it + 
to quantize the system, where again ... denotes the terms out of the round bracket in expressions (20) or (24). Using 
this Hamiltonian in the hermitian form a _1 n^ = n a a _1 n a and also representing n a by —id a , this is equivalent to our 
above treatment in which the last term in (39) is absent. Therefore, one may have some doubts on the validity of the 
main following results due to the effects of the chosen factor ordering. To overcome this problem, we have made some 
calculations based on the above mentioned transformed Hamiltonian and have verified that the general patterns of the 
resulting wave functions follow the behavior shown in following sections. 



9 



This equation does not seem to have analytical solutions. However, we can get some properties of its 
solutions in special regions where there is interest in classical and quantum regimes. First of all, let 
us rewrite this equation in the form 

j-^ - (V 1 + 6iAa 3 ) + [(36 - 6iA) a 2 + 12Aa 4 - 9A 2 a 6 ] j *(o) = 0, (41) 

in which we have used the numerical values C_ = —9/4 and c_ = 3/4 [6]. For large values of 
a, the solution to this equation can easily be obtained in the Wentzel-Kramers-Brillouin (WKB) 
(semiclassical) approximation. In this regime we can neglect the term a^ 1 in equation (41). Then, 
substituting ^(a) = f2(a)e* S(a ) in this equation leads to the modified Hamilton-Jacobi equation 

_ gT) 2 + 6A a 3 ^ + (36a 2 + 12Aa 4 - 9A 2 a 6 ) + Q = 0, (42) 

in which the quantum potential is defined as Q = ^^r- It is well-known that the quantum effects 
are important for small values of the scale factor and in the limit of the large scale factor can be 
neglected. Therefore, in the semiclassical approximation region we can omit the Q term in (42) and 
obtain 

— = 3Aa 3 ± aV36 + 12Aa 2 . (43) 
da 

In the WKB method, the correlation between classical and quantum solutions is given by the relation 
P a = Thus, using the definition of P a in (19), the equation for the classical trajectories becomes 




from which one finds 



(44) 



a(t) = y|sinh l±J±{t-S)\. (Jo) 



which shows that the late time behavior of the classical cosmology (27) is exactly recovered. The 
meaning of this result is that for large values of the scale factor, the effective action corresponding to 
the expanding and contracting universes is very large and the universe can be described classically. 
On the other hand, for small values of the scale factor we cannot neglect the quantum effects, and the 
classical description breaks down. Since the WKB approximation is no longer valid in this regime, 
one should go beyond the semiclassical approximation. In the quantum regime, if we neglect the term 
A 2 a 6 in (41), the two linearly independent solutions to this equation can be expressed in terms of the 
Hermite H v {x) and hypergeometric \F\ (a, b; z) functions, leading to the following general solution: 



*(a) = e 



— ia 2 



/ (l + »)(2 + 3Aa 2 ) \ fl 4 1 i(2 + 3Aa 2 ) 2 " 

C\n_i «. 7= + c 2 1^1 - + —i, -, 



2 



3A l V 2V3A / \4 3A T 6A 



(46) 



At this step we take a quick glance at the question of the boundary conditions on the solutions 
to the Wheeler-DeWitt equation. Note that the minisuperspace of the above model has only one 
degree of freedom denoted by the scale factor a in the range < a < oo. According to [18], its 
nonsingular boundary is the line a = 0, while at the singular boundary this variable is infinite. 
Since the minisuperspace variable is restricted to the abovementioned domain, the minisuperspace 
quantization deals only with wave functions defined on this region. Therefore, to construct the 
quantum version of the model, one should take into account this issue. This is because in such cases, 
one usually has to impose boundary conditions on the allowed wave functions; otherwise the relevant 
operators, especially the Hamiltonian, will not be self-adjoint. The condition for the Hamiltonian 
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operator ~H associated with the classical Hamiltonian function (20) and (24) to be self-adjoint is 

(ipufifo) = (^1,^2) or 

poo poo 

/ ipl(a)Hifo(a)da = / ^ 2 {a)U^l{a)da. (47) 

JO JO 

Following the calculations in [19] and dealing only with square integrable wave functions, this con- 
dition yields a vanishing wave function at the nonsingular boundary of the minisuperspace. Hence, 
we impose the boundary condition on the solutions (46) such that at the nonsingular boundary (at 
a = 0), the wave function vanishes. This makes the Hamiltonian hermitian and self-adjoint and can 
avoid the singularities of the classical theory, i.e. there is zero probability for observing a singularity 
corresponding to a = 0. 2 Therefore, we require 

c 2 H -hm i (v§ 



*(o = 0) = => — = / 3A ivgw ( 48 ) 

11 \^ 4 ^ 3A > 2 > 3A / 

Note that equation (41) is a Schrodinger-like equation for a fictitious particle with zero energy moving 
in the field of the superpotential with the real part U(a) = —(36a 2 + 12Aa 4 ). Usually, in the 
presence of such a potential the minisuperspace can be divided into two regions, U > and U < 0, 
which could be termed the classically forbidden and classically allowed regions, respectively. In the 
classically forbidden region the behavior of the wave function is exponential, while in the classically 
allowed region the wave function behaves oscillatorily. In the quantum tunneling approach [18], 
the wave function is so constructed as to create a universe emerging from nothing by a tunneling 
procedure through a potential barrier in the sense of usual quantum mechanics. Now, in our model, 
the superpotential is always negative, which means that there is no possibility of tunneling anymore, 
since a zero energy system is always above the superpotential. In such a case, tunneling is no 
longer required as classical evolution is possible. As a consequence the wave function always exhibits 
oscillatory behavior. In figure 2, we have plotted the square of the wave functions for typical values 
of the parameters. It is seen from this figure that the wave function has a well-defined behavior near 
a = and describes a universe emerging out of nothing without any tunneling. (See [22], in which such 
wave functions also appeared in the case study of the probability of quantum creation of compact, 
flat, and open de Sitter universes.) On the other hand, the emergence of several peaks in the wave 
function may be interpreted as a representation of different quantum states that may communicate 
with each other through tunneling. This means that there are different possible universes (states) 
from which the present universe could have evolved and tunneled in the past, from one universe 
(state) to another. 

4.2 Perfect fluid quantum cosmology 

In this case, the Wheeler-DeWitt equation can be constructed by means of the Hamiltonian (24). 
With the same approximations as we used in the previous subsection, we obtain 




oa 



(36 - 6iA) a 2 + 12Aa 4 ] - ia 1 ^-^ 1 *(a, T) = 0. (49) 



We separate the variables in this equation as 

tf(a,T) = e iJ52 >(a), (50) 

2 Such a boundary condition is also suggested by DeWitt in the form ^[t/' 3 '] = [20], where C/' 3 ' denotes all three- 
geometries which may play the roll of barriers, for instance singular three-geometries. As is argued in [20], with this 
boundary condition some kinds of classical singularities can be removed and a unique solution to the Wheeler-DeWitt 
equation may be obtained. Although in the presence of more fundamental proposals of the boundary condition in 
quantum cosmology (for example, Vilenkin's tunneling or Hawking's no boundary proposals), it is not clear that the 
above mentioned boundary condition is true, there are some evidences in quantum gravity models in which suitable 
wave packets obey such kind of boundary condition, see [21]. 
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Fi gure 2: The square of the wave function for the quantum universe. We take the numerical value A = 1.5. 



leading to 



|^ - (a" 1 + 6iAa 3 ) + [(36 - 6»A) a 2 + 12Aa 4 + Ea 1 ^] j 

The solutions of the above differential equation may be written in the form 

'(1 + i)(2 + 3Aa 2 Y 



V>(a) = 0. 



(51) 



1> E (a) = e~ ia ' Cl H_ 
for oj = —1/3 and 
tp E (a) = 



1 32+ E ■ 



2V3A 



„ 1 32 + E 1 i(2 + 3Aa 2 ) 2 
+ C2li?1 4 + ^iA- Z '2 ; 6A 



(52) 



-i(l + T fr)a 2 



C2 1-^1 



Cjil x 1152A 2 -24EA-fi 2 
2 432A^ 



'(1 + i) + 6A(2 + 3Ao 2 )] 



12AV3A 



+ 



' 1 + 1152A 2 - 24EA — E 2 . lj[E + 6A(2 + 3Aa 



2M2\1 



^4 ' 864A 3 '2' 216A 3 

for oj = — 1. Now the eigenfunctions of the Wheeler-DeWitt equation can be written as 

* E (a,T) = e iET iP E {a). 



(53) 



(54) 



We may now write the general solution to the Wheeler-DeWitt equation as a superposition of its 
eigenfunctions; that is, 



*(a,T) 



A(E)V E (a,T)dE, 



(55) 



where A(E) is a suitable weight function to construct the wave packets. The above relations seem 
to be too complicated to extract an analytical expression for the wave function. Therefore, in the 
following (for the case u = —1/3), we present an approximate analytic method which is valid for very 
small values of scale factor, i.e., in the range that we expect the quantum effects to be important. 
In this regime if we keep only the a -1 and o 2 terms in the second and third terms of (51), the 
solutions to this equation can be viewed as a superposition of the functions sin ^ v / 36+g- _faA Q 2^ an( j 
cos ^ v / 36+g- _faA Q 2^ jf we i m p 0ge the boundary condition ip(a = 0) = on these solutions, we are led 
to the following eigenfunctions: 

y E (a,T) = e iET sin 

Now, by using the equality 

roo 

e~ lx sin Jmxdx 



V36 + E - 6«A 



(56) 



2^/2' 



-(m/47) 



(57) 
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we can evaluate the integral over E in (55), and the simple analytical expression for this integral is 
found if we choose the function A(E) to be a quasi-Gaussian weight factor A{£ ) = e~ l£ (7 is an 
arbitrary positive constant and £ = 36 + E — 6iA), which results in 



*(a, T) = e -6(A+6*)T 00 e - l£e i£T sin ( V^ a 2 j d£ _ (5g) 



Using the relation (57) yields the following expression for the wave function 

• (o,r, ^«_^_ exp (-^) , ,59) 

where M is a numerical factor. Now, having this expression for the wave function of the universe, we 
are going to obtain the predictions for the behavior of the dynamical variables in the corresponding 
cosmological model. To do this, one may calculate the time dependence of the expectation value of 
a dynamical variable q as 

< VD r lo|\E' > 

< q > (T) = ^¥^- (60) 

Following this approach, we may write the expectation value for the scale factor as 

Jg°**{a,T)*(a,T)da ' ( j 



which yields 



This relation may be interpreted as the quantum counterpart of the classical solutions (34). However, 
in spite of the classical solutions, for the wave function (59), the expectation value (62) of a never 
vanishes, showing that these states are nonsingular. Indeed, in (62) T varies from —00 to +00, and 
any To is just a specific moment without any particular physical meaning like big-bang singularity. 
The above result may be written in terms of the cosmic time r. By the definition dr = a~ 1 (T)dT, 
we obtain the quantum version of the relations (35) as 



< a > (r) = i (e^ T + 7 2 e - V^ T ^ . (63) 

In figure 3, we have plotted the classical scale factors (34) and (35) and their quantum counterparts 
(62) and (63). As is clear from this figure, for a perfect fluid with oj = —1/3, the corresponding classical 
cosmology admits two separate solutions which are disconnected from each other by a classically 
forbidden region. One of these solutions represents a contracting universe ending in a singularity 
while another describes an expanding universe which begins its evolution with a big-bang singularity. 
On the other hand, the evolution of the scale factor based on the quantum-mechanical considerations 
shows a bouncing behavior in which the universe bounces from a contraction epoch to a reexpansion 
era. Indeed, the classically forbidden region is where the quantum bounce has occurred. We see that 
in the late time of cosmic evolution in which the quantum effects are negligible, these two behaviors 
coincide with each other. This means that the quantum structure which we have constructed has a 
good correlation with its classical counterpart. 



5 Bohmian trajectories 

In the previous sections, we saw how the classical singular behavior of the universe was replaced with a 
bouncing one in a quantum picture. Now, a natural question may arise: Why will the bounce occur? 



13 



Fi gure 3: Left: The figure shows qualitative behavior of the classical scale factor (34) (solid lines: the left branch 
for a^\t) and the right branch for a fII \t)) and the expectation value of the scale factor (62) (dashed line). Right: 
The same figure in terms of cosmic time. The left and right branches of the solid lines represent a^(r) and a' J/ '(r), 
respectively, in (35) while the dashed line represents the expectation value (63). 



Clearly, it is due to the quantum mechanical effects which show themselves when the size of the 
universe tends to very small values. However, we would like to know whether the massive correction 
to the underlying gravity theory has any contribution to this phenomenon. To deal with this question, 
let us return to the wave function (59) and write it in the polar form ^(a, T) = 0(a, T)e lS ( a ' T \ where 
U(a, T) and S(a, T) are real functions, which simple algebra gives as 



7a 2 



8( 7 2 + T 2 ) 



(64) 



3 T cl^ 

S(a, T) = -36T + - arctan - - g( 2+r2) - (65) 

According to the Bohm-de Broglie interpretation of quantum mechanics [23] and also its usage in 
quantum cosmology [24], upon using this form of the wave function in the corresponding wave equa- 
tion, we arrive at the modified Hamilton-Jacobi equation as 

n { quPi = w) + Q = ' ( 66 ) 

where Pi are the momentum conjugate to the dynamical variables qi and Q is the quantum potential. 
With beginning of the wave equation (49), for which we have used the same approximations as in the 
previous section, the above mentioned procedure gives the quantum potential as 

1_^_J_^ 

^ nda 2 aQda V ' 

On the other hand, the Bohmian equations of motion can be obtained by P a = where by means 
of the relation (19) reads 

T 

-6aa + 3Aa 2 = (68) 
The solution to this equation denotes the Bohmian representation of the scale factor; that is 



a(t) = ^Jce AT + ^e AT "^ A [e 2 ^ A Ei(l; -AT - i 7 A) + Ei(l; — AT + irfA)], (69) 
where c is an integration constant and Ei(6; z) is the exponential integral function defined by 

/oo 
e~ kz k- b dk. (70) 

The bouncing behavior of the scale factor is again its main property near the classical singularities as 
we have shown in figure 4. To achieve an expression for the quantum potential in terms of the scale 
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a(Bolimian) 




a 



Figure 4: Left: The Bohmian trajectory of the scale factor. Right: The horizontal line represents a typical energy 
level. The solid curve is the quantum potential for A / 0, while the dashed curve denotes the quantum potential when 
A = 0. 



factor, we note that all of our above calculations are in the vicinity of T ~ 0, where the scale factor 
is small. In this regime, a numerical analysis shows that the Bohmian scale factor (69) behaves as 
a(T) ~ (7 2 + T 2 ) 1 / 2 , in agreement with the expectation value (62). Thus, substituting in (64), we 
get the quantum potential from (67) as 

Q(«) = \ 

In figure 4 we also have plotted the qualitative behavior of the quantum potential versus the 
scale factor. As this figure shows, this potential goes to zero for the large values of the scale factor. 
This behavior is expected, since in this regime the quantum effects can be neglected and the universe 
evolves classically. On the other hand, for the small values of the scale factor the potential takes a 
large magnitude and the quantum mechanical considerations come into the scenario. This is where the 
quantum potential can produce a huge repulsive force which may be interpreted as being responsible 
of the avoidance of singularity. In figure 4 the horizontal line represents a constant energy level which 
in intersecting with the potential curves gives the turning points at which the bounce will occur. The 
solid curve in this figure is plotted in the case of A ^ 0; i.e., for the massive theory, while the dashed 
curve is for A = 0; i.e., for when the massive corrections are absent. It is seen that, although the mass 
term A is not the only reason for the bouncing behavior in the vicinity of the classical singularity, 
it may shift the bouncing point into the smaller values of the scale factor. This means that if we 
consider the bouncing point as the minimum size of the universe (which is suggested by quantum 
cosmology), then the massive version of the underlying gravity theory predicts a smaller value for 
this minimal size in comparison with the usual Einstein-Hilbert model. These facts and also other 
considerable possibilities such as quantum tunneling between different classically allowed regimes (as 
can be seen from figure 4) through the potential barrier support the idea that the massive corrections 
to the classical cosmology are some signals from quantum gravity. 



+ 



8A 



a*- 7 2 a 2 ' (a 2 - 7 2 ) 3 / 2 



+ 



-1 + 48AV 



(71) 



6 Conclusions 

In this paper we have applied the recently proposed nonlinear massive theory of gravity to an open 
FRW cosmological setting. Although the absence of homogeneous and isotropic solutions is one of 
the main challenges related to this kind of gravitational theory, we moved along the lines of [6, 7], in 
which the existence of open FRW cosmologies is investigated. By using the constraint corresponding 
to the Stiickelberg scalars, we reduced the number of degrees of freedom, according to which the 
total Hamiltonian of the model is deduced. We then presented in detail, the classical cosmological 
solutions either for the empty universe or in the case where the universe is filled by a perfect fluid 
(in its Schutz representation) with the equation of state parameter cj = —1/3, —1. We saw that in 
both of these cases, the solutions consist of a contraction universe which finalizes its evolution in 
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a singular point and an expanding universe which begins its dynamic with a big-bang singularity. 
These two branches of solutions are disconnected from each other by a classically forbidden region. 
Also, the common feature of the vacuum and matter classical solutions is that the mass term plays 
a role which resembles the role of cosmological constant in the usual de Sitter universe. In this sense 
we may relate the massive corrections of GR to the problem of dark energy. 

In another part of the paper, we dealt with the quantization of the model described above via 
the method of canonical quantization. For an empty universe, we have shown that by applying the 
WKB approximation on the Wheeler-DeWitt equation, one can recover the late time behavior of the 
classical solutions. For the early universe, we obtained oscillatory quantum states free of classical 
singularities by which two branches of classical solutions may communicate with each other. In the 
presence of matter, we focused our attention on the approximate analytical solutions to the Wheeler- 
DeWitt equation in the domain of small scale factor, i.e. in the region which the quantum cosmology is 
expected to be dominant. Using Schutzs representation for the perfect fluid, under a particular gauge 
choice, we led to the identification of a time parameter which allowed us to study the time evolution 
of the resulting wave function. Investigation of the expectation value of the scale factor shows a 
bouncing behavior near the classical singularity. In addition to singularity avoidance, the appearance 
of bounce in the quantum model is also interesting in its nature due to prediction of a minimal size for 
the corresponding universe. We know the idea of existence of a minimal length in nature is supported 
by almost all candidates of quantum gravity. Finally, we repeated the quantum calculations by means 
of the Bohmian approach to quantum mechanics. The analysis of the quantum potential shows the 
importance of the mass term in the action of the model. Indeed, we have shown that in the presence 
of the massive graviton, the quantum potential changes its behavior from an infinite barrier to a finite 
one, and hence the minimal size of the universe, from which the bounce occurs, will be shifted to 
the smaller values. Also, the massive theory of quantum cosmology exhibits some other possibilities; 
for example, tunneling between different classically allowed regions, for cosmic evolution in the early 
universe epoch. 
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